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ABSTRACT 

New convenient stability criteria are provided in this paper for a 
large class of finite difference approximations to initial-boundary 
value problems associated with the h}^erbolic system 
-t- Bm -f f in the quarter plane x > 0, f > 0. Using the new 
criteria, stability is easily established for numerous combinations of 
well known basic schemes and boundary conditions, thus generaliz- 
ing many special cases studied in recent literature. 
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0. Introduction. 


In this paper we extend the results of [3] to obtain convenient, more versa- 
tile, sufficient stability criteria for a wide class of difference approximations for 
initial-boundary value problems associated with the hyperbolic system 
u< = Awa + + f in the quarter plane x > 0, f > 0. Our difference approxima- 

tion consists of a general difference scheme - explicit or implicit, dissipative or 
not, two-level or multi-level -- and boundary conditions of the type discussed in 
[3]. 

We restrict attention to the case where the outflow part of our boundary 
conditions is translatory; i.e., determined at all boundary points by the same 
coefficients. Such boundary conditions are commonly used in practice; and in 
particular, when the boundary consists of a single point, the boundary condi- 
tions are treinslatory by definition. 

Throughout the paper we assume that the basic scheme is stable for the 
pure Cauchy problem, and that the assumptions that guarantee the validity of 
the stability theory of Gustafsson, IO*eiss and Sundstrom [5] hold in our case. 
With this in mind, we raise the question of stability for the entire approximation 
in the sense of Definition 3.3 in [5]. 

We begin our stability analysis in Section 1 by recalling Theorem 2.1 of [3], 
which implies that the entire approximation is stable if and only if the scalar 
outflow components of its principal peirt are stable. Thus, our global stability 
question is reduced to that of a scalar, homogeneous, outflow problem which is 
the mean subject of this paper. 

Our stability criteria for the reduced problem -- stated in Section 2 and pro- 
ven in Section 3 — depend both on the basic scheme and the boundary condi- 
tions, but very little on the intricate interaction between the two. Moreover, we 
show in Remark 3.1 that our old results in [3] easily follow from the present 
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ones. Thus, our new criteria provide in many cases a convenient alternative to 
the well known Gustafsson-Kreiss-Sundstrom criterion in [5], 

In the second part of Section 2 we use the new stability criteria to reestab- 
lish all the main examples in our previous paper [S], We show, for instance, that 
if the basic scheme is arbitrary (dissipative or not) and the boundary conditions 
are generated by either the explicit or implicit right-sided Euler schemes, then 
overall stability is assured. For dissipative basic schemes we prove stability if 
the boundary conditions are determined by either oblique extrapolation, by the 
Box-Scheme, or by the right-sided weighted Euler scheme. These and other 
examples incorporate most of the cases discussed in recent literature [l]- 
[3],[5].[6],[8]-[l0].[l2].[l3].[l5].[l6]. 

Section 2 contains some new examples as well. Among these we find that if 
the basic scheme is arbitrary and two-level, then horizontal extrapolation at the 
boundary maintains overall stability. Other stable cases occur when the basic 
scheme is given by either the backward (implicit) Euler scheme or by the 
Crank-Nicolson scheme and the boundary conditions are determined by oblique 
extrapolation. Such examples, where neither the basic scheme nor the boun- 
dary conditions are necessarily dissipative, could not have been handled by our 
previous results in [3]. 

As in [3], we point out that there is no difficulty to extend our stability cri- 
teria to problems with two boundaries. This is so since if the corresponding left 
and right quarter-plane problems are stable, then by Theorem 5.4 of [5], the ori- 
ginal problem is stable as well. 
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1. The Difference i^^roximation and the Reduced Problem. 

Consider the first order h)^erbolic system of partial differential equations 


9u(x, t)/dt = /49u(x, t)/dx + 5u(x, t) + f(x, t), x > 0. f > 0. (1.1a) 


where u(x, t) = (u^*)(x, t), . . . t))' is the unknown vector (prime denot- 
ing the transpose), t(x, t) = t) /^”^(x, t))' is a given vector, and A 

and B are fixed nxn matrices such that A is Hermitian and nonsingular. Without 
restriction we may assume that A is diagonal of the form 


A = 


0 


0 

• 


A^ < 0 . A^^ > 0 , 


( 1 . 2 ) 


where A^ and A^^ are of orders Ixi and (n — Z) x (n — 1), respectively. 

The solution of (1.1) is determined uniquely if we prescribe initial values 


u(x, 0) = u°(x), X > 0, (1-lb) 

cind boundary conditions 

u^(0, 0 = S'u^^(0. Z) + g(0, (1-lc) 

where 5 is a fixed I x(n, — 1) matrix, g(Z ) is a given Z-vector, and 

u^ = (u^') u^^ = uf"))' (1-3) 

is a partition of u into its inflow and outflow components, corresponding to the 
partition of ^4 in (1.2). 

In order to solve the initial-boundary value problem (l.l) by a finite 
difference approximation we introduce, as usual, a mesh size Ax > 0, AZ > 0, 
such that X = AZ/Ax = constant: and using the notation Vj,(Z) = v(i/Ax, Z), we 
approximate (1.1a) by a general, consistent, two-sided, solvable, multi-level 


basic scheme of the form 
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Q-iVy{t + Ai) = S Qc^vit - ctAO + 

< 7=0 

(1.4a) 

Qa= ^ AjoE^ , £’V^ = V„+ 1 , CT = -1 s, i/ = r,r + l 

j=-r 

Here, the tixti matrices Aj^j are pol}niomials in A and ^tB, and the n-vectors 
bi/(^) depend smoothly on f(x, t) and its derivatives. 

The equations in (1.4a) have a unique solution if we provide initial values 

Vj,(aAO = v^(aAO. a = 0 s. i/ = 0, 1, 2 (1.4b) 

where in addition we must specify, at each time step t = //Af , = s, s + 1, . . . , 

boundary values v^(f + Af). i/ = 0, . . . ,r - 1. As in [3], these boundary values 
will be determined by two sets of boundary conditions, the first of which is 
obtained by taking the last n — I components of general boundary conditions of 
the form 

T-iVy{t + AO = ^ - aM) + A^d^(^). 

< 7=0 


Ta= t Cj„E^, a=-l g, u = 0 r - 1, 

i=o 


where the matrices Ov ^ire pol5momials in A eind. MB , the Cj(-i) nonsingular, 
and the n-vectors dj,(<) are functions of f{x, t) and its derivatives. 

If we put 





in accordance with the partition of A and u in (1.2), (1.3), then this set of condi- 
tions takes the explicit form 
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r?/ vlit + AO + + AO 

= ^ vl{t - aAO + vl^it - ctAoI + Mdl^t). (1.4c) 

<7=01 J 

Ti^‘" = 2 . a. = I, II, v = Q r - 1. 

j=0 

We call such boundary conditions translatory since they are determined at all 
boundary points by the same coefficients. 

For the second set of boundary conditions we use the analytic condition 

Vo(^ + Af) = SvQ{t + AO + g{t + A^), (l-4d) 

together with r - 1 additional conditions (not necessarily translatory) of the 
form 

+ AO = 2 [ -D// v/(i + Af) + DjJ^ v/(i + AO] + Aie^(0. 
i=i 

(l.4e) 

u = I r - 1. 

where the matrices DjJ and dJJ^ —of orders Ix-l and lx{n — 1) respectively — 
are pol 3 Tiomials in the blocks and a, ^ = I, II, of the matching parti- 

tions 



so that the DjJ vanish whenever B does; and the i-vectors Bl{t) are functions of 
f(x, 0. g(0> their derivatives. 

It has been shown in [3] that equations (l.4c, d, e) can be solved uniquely 
for the required boundary values vjt + AO. = 0 r - 1, in terms of neigh- 

boring values of v, at least for sufficiently small Af; and that boundary condi- 
tions of the form (1.4e) can be constructed to any degree of accuracy. A con- 
crete example of second order accurate boundary conditions of the form (1.4e), 
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for the special case B = f = 0. was given in [2], 

The difference approximation is completely defined now by (1.4); so assum- 
ing that the basic scheme is stable for the pure Cauchy problem (— ~ < i/ < “), 
we may ask whether the entire approximation is stable. More precisely, we 
make from now on the same assumptions about approximation (1.4) as in [5], so 
that the stability theory of Gustafsson, Kreiss and Sundstrom holds, and raise 
the above stability question in the sense of Definition 3.3 of [5]. 

As in [3], the first step will be to reduce our stability question to that of a 
scalar, outflow approximation with homogeneous translatory boundary condi- 
tions. This reduction is obtained by applying (1.4a, b, c) to the scalar outflow 
problem 

Ut = aUj.. a = constant >0, x >: 0, t > 0 (1-5) 

(which requires no analytic boundary conditions). In other words, we set 

= B = f = 0. = a = constant > 0; 


so that (l.la) yields (1.5); and (1.4a, b, c) reduces to a self-contained, scalar, 
homogeneous approximation which consists of the basic scheme 


+ Af) = £ Q„Vy{t — aAf). v = r , r + I, , 

( 7=0 

(1.6a) 

Qq= a = -1, . . . ,s; 

i=-r 


with initial values 


v^,{aLt) -=v^{aht), ct = 0 s, = 0, 1, 2,... ; (1.6b) 


and translatory outflow boundary conditions 
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+ td) = TaVv{t - ctAO. II = 0 r — 1, 

c=0 

(1.6c) 

m 

T„= I, Cj„E> , Co(-i) 0, (7 = -1 q. 

j=0 

where the scalars Ojg and Cjg are polynomials in a, and the basic scheme (1.6a) 
is consistent with (1.5). 

We are now ready to use our main result in Section 2 of [3], which we refor- 
mulate as follows; 

THEOREM 1.1 (Theorem 2.1. [3]). Approximation (1.4) is stable if and only 
if its reduced form in. (1.6) is stable. 

In other words, we have that approximation (1.4) is stable if and only if the 
scalar outflow components of its principal part are stable. 

The above reduction theorem implies that from now on we may restrict our 
stability study to the scalar outflow approximation (1.6). Hence, we conclude 
this section by stating the following five assumptions which will hold throughout 
the paper, and guarantee the vedidity of the Gustafsson-Kreiss-Sundstrom theory 
[5] for approximation (1.6). 

Assumption 1.1 (Assumption 3.1, [5]; Assumption 2.1, [3]). Approximation 
(1.6) is boundedly solvable; i.e., there exists a constant K > Q such that for each 
2 /£ lz{^) there is a unique solution rue Izibx) to the equations 

Q-iw^, = yv. v = r,r + 1 

T-lW^, = Vk . 1 / = 0 r - 1, 

with ||ru|| ^ K\\y\\ ; where 9-i smd T-\ are defined in (1.6a, c), and lz{C^) is the 
space of all grid functions w = (ro^J“=o with ||ru||® = Ax 2"=o < “• 

Assumption 1.2 (Assumption 5.1, [5]; Assumption 2.2, [3]). The basic 
scheme is stable for the pure Cauchy problem -“ < i/ < ~. That is, if we define 
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the basic characteristic function \>y 

P{z,ic)= Oj{z)ic^ (1.7) 

j=~r 

where 

aj{z) s - i j = -r. ... ,p. (l.B) 

( 7=0 

then we have: 

(i) The von Neumann condition; i.e.. the solutions z(/c) of the basic charac- 
teristic equation 

P{z,k) = 0 (1.9) 

satisfy 

|z(/c)| < 1 for all ic'with |/c| = 1. 

(ii) If |/c| = 1 , and if z(k) is a root of (1.9) with \z{k) \ = 1, then z(/c) is a 
simple root of (1.9). 

Assumption 1.3 (see Assumption 5.4 together with Definition 10.1, [5]; com- 
pare Osher [10]). The basic scheme (1.6a) belongs to the family of schemes for 
which the Gustafsson-IQ’eiss-Sundstrom theory in [5] holds. This family contains 
in particular the following two classes: 

(i) Dissipative basic schemes; i.e., schemes for which the roots z(/c) of 
(1.9) satisfy 

|z(/c)|<l for all KUjith\K\ = 1, K i, (1.10) 

(ii) Unitary basic schemes (also known as strictly nondissipative 
schemes) where the roots of (1.9) satisfy 

|z(/c)|=l forall\K\=l. (1.11) 

Obviously, if the basic scheme belongs to any of these two classes, then it 
satisfies the von Neumann condition in Assumption 1.2 (i). 
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Assumption 1.4 (Assumption 5.5. [5]; Assumption 2.4, [3]). 

a_r(z), Op{z) ^ 0 /or all\z\>. 1. 

Assumption 1.5. We assume that 

2 |c,-( 2 )| 0 for oii |z I & 1 

,=o 

where, in analogy with (1.8), 

^/(^) ^ ~ ^ ^ ^ ' j ~ 0, ,, , ,m. (1.12) 

c=o 

Assumption 1.5 is necessary for stability, as shown in Remark 3.4 below. 
This last assumption — which should have been included in [3] as well — is easily 
verified for cdl practical boundary conditions. 

2. Statement of Results and Examples. 

In order to state our main stability criteria we define, in complete analogy 
with (1.7), the boundary characteristic function 

J V ^ y V • 

R(z, k) = 2 cAz)k^ 

3=0 

where the c^ (z) are given in (1.12). Defining the function 

n(z, /c) s |P(z, k)\ + \R{z, /c)| 
we shall prove in Section 3: 

THEOREM 2.1 (1st Main Theorem). Approximation (1.6) is stable if 

n(z, k) 0 for all |z | & 1, 0 < |k| ^ 1, (z, k) ^ (1, 1). (2.1) 

Next, let us divide the (z, /c) domain in (2.1) into three disjoint parts and 
restate Theorem 2.1 as follows: 
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THEOREM 2.r. (1st Main Theorem revisited) ^aproximation (1.6) is stable if 
n(z, k) 7^ 0 for all |z | > 1, |k| = 1, /c 1; (2.2a) 

n(z. /c = 1) 0 for ail |z | ^ 1, z 1; (2.2b) 

n(z, «) 0 for ail |z I > 1, 0 < |«:| < 1. (2.2c) 

The advantage of this setting is explained by Theorem 2.2 in which we pro- 
vide useful suflicient conditions for each of the three inequalities in (2.2) to hold. 
Before stating this theorem, however, we need the following definitions: 

Definition 2.1, The boundary conditions (1.6c) are said to be dissipative if 
the roots z (k) of the boundary characteristic equation 

R{z,ic) = Q (2.3) 

satisfy 

| z ( k )|<1 for ail\K\ < 1, k 7^ 1. 

Definition 2.2. We say that the boundary conditions (1.6c) satisfy the von 
Neumann condition if the roots z(/c) of (2.3) satisfy 

I z ( k) I < 1 for ail\ic\ = 1 . 

Definition 2.3. The boundary conditions (1.6c) are called boundediy solv- 
able if there exists a constant A" > 0 so that for each y £ f 2 (^) there is a unique 
solution iw£ fa(^) to 

T-iWy = yy, v= 0, 1,2 (2.4) 

withll-ujll^ A'ilyll. 

Evidently, these definitions are analogous to those made for the basic 
scheme in Assumptions 1.3(i), 1.2(i) and 1.1, respectively: and again, dissipa- 
tivity implies the von Neumann condition. 
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Having these definitions we shall prove: 

THEOREM 2,2 (2nd Main Theorem). 

(i) If either the basic scheme or the ho'UJidary conditions are dissipative, 
then (2.2a) holds, 

(ii) Inequality (2.2b) is satisfied if any one of the following holds: 

(a) The basic scheme is two-level. 

(b) The basic scheme is three-level and 

n(z = -1. /c= 1) 5^ 0. (2.5) 

(c) The boundary conditions are two-level and at least of zero order of 
accuracy. 

(d) The boundary conditions are three-level, at least zero order accu- 
rate, and (2.5) is satisfied. 

(iii) If the boundary condilioT^s satisfy the von Neumann condition and are 
boundedly solvable, then (2.2c) holds. 

The proofs of Theorems 2.1 and 2.2 will be given in Section 3. 

We turn now to examples. Before doing so, however, let us recall one more 
result which provides a useful sufficient condition for the boundary conditions to 
be boundedly solvable. 

LEMMA 2.1 (Lemma 3.2, [3]). 

(i) The boundai^/ conditions (1.6c) are boundedly solvable if and only if 

T^i{k) = ^ ^ 0 for aZZ 0 < |a:| < 1. (2.6) 

y=o 

(ii) In particular, explicit boundary conditioTis are always boundedly solv- 
able. 

This result is associated with important observations on solvability by Osher 


[ 11 ]. 
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We remark that part (i) of Lemma 3.2 in [3] merely states that if (2.6) holds, 
then the boundary conditions are boundedly solvable. Following the proof of 
Lemma 3.2 in [3], however, one realizes without difficulty that in fact, (2.6) is 
both sufficient and necessary for the bounded solvability of the boundary condi- 
tions (1.6c), as stated above. 

Example 2.1 ([3, Example 3.5]; compare the special cases in [5, Theorems 
6.1, 6.2 and 6.3], [13, (3.2)], and [2, Example 2]). Let the basic scheme (1.6a) be 
arbitrary (dissipative or not), and let the boundary conditions (1.6c) be gen- 
erated by the right-sided, first order accurate, explicit Euler scheme: 

v^,{t + At) = + Xa[vt,^.i{t) -v^,{t)], 0 < Xo < 1, u=0 r - 1. 

These two-level boundary conditions eire known to be dissipative (see Example 
3.5, [3]), hence Theorem 2.2(i)(iic) implies (2.2a, b). Since the boundary condi- 
tions are explicit, then by Lemma 2.1 they are solvable: and since they are dissi- 
pative, the von Neumann condition is satisfied. Consequently, Theorem 2.2(iii) 
3 delds (2.2c), and stability is established by Theorem 2. 1'. 

It is a trivial matter to verify that for this, as well as for all the following 
excimples. Assumption 1.5 is fulfilled. 

Example 2.2 ([3, Example 3.6]; compare the special cases [13, (3.3)], and [2, 
Example 3]). Take an arbitrary basic scheme, and determine the boundary con- 
dition by the right-sided, first order accurate, implicit Euler scheme: 

Vi,{t + At) - Xa[v,/+i(f + At) — Vf,{t + Af)] = ^^(f), Xa > 0, v = 0 r — 1. 

As in the previous example, the boundary conditions are two-level and dissipa- 
tive (e.g. Example 3.6, [3]), so Theorem 2.2(i)(iic) gives (2.2a, b). The dissipa- 
tivity of the boundary conditions also implies the von Neumann condition. And it 
is trivially verified (Example 3.6 [3] again) that 7’_i(«:) 0 for \ic\ ^ 1 (where 

T-i{k) is defined in (2.6)), so that Lemma 2.1(i) implies solvability. Thus, (2.2c) 
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follows from Theorem 2.2(iii); and Theorem 2.1' assures stability. 

ExamplB 2.3 (Compare the special cases in [B, Theorem 6], [10. Section 
XXIII], [5, Theorems 6.1 and 6.3], [1], [6. Theorem 2.1], and [3. Example 3.1]). 
Take an arbitrary two-level basic scheme, and define the boundary conditions by 
horizontal extrapolation of order u - 1: 

+ AO = 2 + AO. = 0 ’’ - 1- (2.7) 

i=i J 

We have 

R{z. k) = i?(o = 1 - s = (1 - o“. 

so R{ic) 9^ 0 for K P 1, which directly gives (2.2a, c). Moreover, since the basic 
scheme is two-level. Theorem 2.2(iia) implies (2.2b); and Theorem 2.1' again 
proves stability. 

It is interesting to note that the above result may fail, both for dissipative 
and nondissipative basic schemes, if the basic scheme is of more than two time- 
levels. A nondissipative counterexample was given by Gustafsson et al. in 
Theorem 6.2 of [5] (see also [6]) who showed that the unitary. 3-level Leap-Frog 
scheme 

v^,{t + A^) = Vt,{t — AO + Xa[v^, + i(0 “ ■*^i.-i(0] 1/ = 1. 2, 3 (2.8) 

provides an unstable approximation in combination with the linear boundary 
extrapolation ((2.7) with « = 2 ): 

VQ{t + At) = 2vi{t + AO “ + AO- 

We proved instability ([3], (3.6)) for the case where the basic scheme is the 3- 
level, 5-point, dissipative version of (2.8): 
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v^{t + AO = [/ - ^{E - If {I - E-^f]v^{t -Lt)+\a{E- E~'^)v^{t), 

0<£<1, 0<Xa<l-£, i/ = 2, 3 

and the boundary conditions are given by (2.7) with i/ = 0, 1. 

Example 2.4 ([3. Example 3.2]; compare the special cases in [2, Example 1] 
and [6, Theorem 2.2]). Let the basic scheme be dissipative, and determine the 
boundary conditions by oblique extrapolation of order a - 1: 

v,{t + AO = - (;• - l)Ai ], 1 / = 0 r - 1. (2.9) 

Since the basic scheme is dissipative, then Theorem 2.2(i) implies (2.2a). Furth- 
ermore, the boundary characteristic function for (2.9) is 

i?(z, 0 = 1 - 2 (?)(-l)^'^'2-V = (1 - Z-‘«)“ 

J=1 J 

hence 

n(z, /c) ^ 01 5^ 0 , z 7i K. ( 2 . 10 ) 

This yields (2.2b, c); and Theorem 2.1’ implies stability. 

Example 2.5 ([3, Example 3.3]; compare the special cases [5, Theorem 6.1], 
[12, (3.4)], and [2, Example 4]). Take ciny dissipative basic scheme, and let the 
boundary conditions be generated by the second order accurate Box-Scheme: 

vj^t + AO + ■i'v+i(^ + Ai) — Aa[^»y+l(^ + AO ~ 'Wi/(^ + Ai)] 

= TJ^(0 +^'u^-l(0 + >'“bv+i(0 -^^ 1 ^( 0 ]. V = 0 r - 1. 

As in the previous example, dissipativity implies (2.2a); and Theorem 2.2(iic) 
implies (2.2b). Next, we recall Example 3.3 of [3] where it was shown without 
difficulty that the roots z{k) of the boundary characteristic function satisfy 
|z(OI = 1 for |/c| = 1. and that r_i(/c) 0 for |/c| ^ 1. Thus, the boundary con- 
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ditions satisfy the von Neumann condition; and Lemma 2.1(i) implies solvability. 
With this, Theorem 2,2(iii) gives (2.2c), and by Theorem 2.1' stability follows. 

Example 2.6 ([3, Example 3.4]; compare the special case in [13. (6.2b)]). 
Take any dissipative basic scheme, and define the boundeiry conditions by the 
right-sided, 3-level, weighted Euler scheme: 

Vy{t + AO = “ ^0 + Xa[27;„^-l(^) - ^'v(^ + AO - vj^t - AO]. 

( 2 . 11 ) 

0 < Xa ^ 1, 1 / = 0 r — 1. 

As in the two previous examples. Theorem 2.2(i) implies (2.2a). Further, we have 
R{z, /c) = 1 - z"® - Xa(2Kz“^ - 1 - 
so 

n(z = -1, c = 1) ^ |^(z = -1, K = 1) I = 4Xa > 0, 

and Theorem 2.2(iid) yields (2.2b). Next, as in Example 3.4 of [3], we find that 
the roots of i?(z, k) satisfy 

z(«: = e'f) = e^f^|^^^, 6(0 = V(Xa)« + e-^f[l - (Xa)“]; (2.12) 

and since 0 < Xa < 1, then 1 6 (f) | ^ 1 for 1 0 ^ tt; therefore 

< 1 . 10 < 77 , 

so the von Ne umann condition holds. By Lemma 2.1(ii) we also have solvability. 
Hence Theorem 2.2(iii) yields (2.2c); and Theorem 2.1' assures stability. 

Example 2.7. Let us keep the boundary conditions in (2.11), and extend the 
result of Example 2.6 to nondissipative basic schemes whose characteristic func- 
tions satisfy 

P(z = -1, ic = -1) 7^ 0. (2.13) 

Evidently, we obtain (2.2b, c) precisely as for the previous example. In order to 
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prove (S.2a) we set p= (Xa)® and observe that b (^) in (2.12) satisfies for 
0 < ll'l < T7, 

|b(OI® =P® + (1 -pf + 2p(l -p) cos 2^ <p® + (1 -pf + 2p(l -p) = 1. 
Thus, by (2.12), 

jz(/c = 1 < I < 1. 0<|^|<7 t; 

or in other words 


R{z, /c) 5^ 0. |z I & 1. |c| = 1, /c ±1. 

By (2.12) again, 


z{k= — 1 ) = 


Xtt ± 1 
Xu + 1 


(2.14) 


hence one root is z (/c = — l) = —1 and the other satisfies \z{k = —1) | < 1; so 

R{z . K = -1) 0, |z|^l, z -1. (2.15) 

Collecting (2.13) - (2.15) we find therefore that 

n(z, k ) = |P(z, c)| + |l?(z, k)| 0, |z I > 1. |/c| = 1, K 1. 

Thus, (2.2a) is established, and Theorem 2. 1 ’ implies stability. 

We remark that certain well known nondissipative basic schemes satisfy 
(2.13). This includes, for example, the unitary Crank-Nicolson scheme 

Vt,{t + At) — ^-[7;,^i(f + Af ) - + Af)] 


(2.16) 


“ i^i/-i(0]. = 1. 2, 3.... : 


and the almost-dissipative* backward Euler scheme 


♦We call a basic scheme almost dissipative if (l.lO) holds, except for a finite number of 
K, Ik I = 1. 
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V^,{t + ^t) — ^-[Vv+l(^ + ^0 “ + ^0] - 

(2.17) 

V = r,r + 1 ,.... 

which is included in the family of schemes for which the Gustafsson-Kreiss- 
Sundstrom theory holds. 

If we consider, however, the almost-dissipative Lax-Friedrichs scheme 

Vy{t + M) = |-[■Uw+l(0 + + ^[^v+i(0 -■yi/-i(0]. 

0 < Xa < 1, V = 1,2, 3,-" 

(which is also admitted by the Gustafsson-Kreiss-Sundstrom theory; e.g. Osher 
[10]), we find that P{z = -i. k = -1) = 0, i.e., (2.13) fails. Indeed, it can be 
shown that this scheme together with the boundary conditions in (2.11) provides 
an unstable approximation. 

Example 2.8 ([15. Section 3]; compare [16. Example 6.2]). Take any two- 
level basic scheme, and determine the boundary conditions by the right-sided, 
two-step Euler scheme: 

vjj: + A^) = v„(f — At) + 2Xa[^^^,^.l(^ -Af)— i/„(f — At)], 

0 < Xtt < 1/ 2, 1 / = 0 r - 1. 

Then, 

R{z, /c) = 1 - z~\l + 2Xa(/c - l)], 

so the roots z (ic) of R{z ,k) satisfy 

I z (k = 1^=1 + 2Xa(2\a — l)(l — cos ^) < 1, 0 < | ^ | ^ rr. 

That is, our explicit boundary conditions are dissipative and satisfy the von Neu- 
mann condition. Hence, by Lemma 2.1(ii) and Theorem 2.2 we obtain (2.2), and 
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Theorem 2.r assures stability. 

Example 2,9 (in [5. Theorem 6.3] and [6. Theorem 2.2]). Consider the uni- 
tary, unconditionally stable Crank-Nicolson scheme in (2.16), with oblique extra- 
polation at the boundary: 

VQ{t + ht) = - (;• - 1)A^] (2.18) 

as in (2.9). Since we have (2.10), then in order to comply with the assumptions 
of Theorem 2.1 it remains to show that 


n(z, «) s |P(z, /c)| + |/?(z, k)| 0, z = K, \K\=l,Kyt\\ 

so it suffices to prove that 


Xa 


P{z, tc) = 1- ^(k -k ^) - Z“l[l + ^{k - K~^)] 0 . 


Z = K = e*f, 0 < 1^1 < 7T. 


(2.19) 


Indeed, 


P(z = e'f, K = e<f) = 1 - cos I + ^(cos® ^ - cos ^ - 1) + i sin ^(l - ^cos 0 , 

so observing that the real and imaginary parts do not vanish simultaneously, we 
immediately obtain (2.19), and stability follows. 

Example 2,10 (in [6, Theorem 2.2]). Take the edmost-dissipative, uncondi- 
tionally stable backward Euler scheme in (2. 17) with the oblique boundary extra- 
polation in (2.18). As in the previous example, in order to prove stability via 
Theorem 2. 1 it is sufficient to show that 

P(z, k) S 1 - ^(/C-/C"‘) - 2 "* ^ 0, z = K = , 0 < 1^1 < 7T. (2.20) 

Since 


P{z — e'’^, K - e'^) = 1 — Xa cos f — cos + i sin |, 
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then, comparing the real arid imaginary parts again, we get (S.20), and stability 
is established. 

At the end of Section 3 we make several remarks concerning these exam- 
ples. 


3. Proofs and Remarks. 

In order to prove Theorems 2.1 and 2.2 we begin by considering the basic 
characteristic function P{z, k) in (1.7). By Assumption 1.4, for each z with 
|z I > 1, P(z , k) has r +p roots /c(z). These roots, which play a major role in 
the stability analysis of approximation (1.6), have the following separation pro- 
perty: 

LEMMA 3.1 ([3], Lemma 4.2; compare Lemmas 5.1 and 5.2, [5]). For \z\ > 1, 
the basic characteristic function P{z,k) has precisely r roots k{z) with 
0 < |/c(z) I < 1, p roots with \ k(z) | > 1, and no roots with |k(z) 1 = 1. 

According to this lemma, the roots k(z) of P{z, k) split for |z | > 1 into two 
groups: r inner roots satisf 3 dng \k{z)\ < 1, and p outer roots with |/c(z)| > 1. 
Using Assumption 1.4 and a continuity argument, we see that these groups of 
inner and outer roots remain well defined for | z | ^ 1 as well, where milder ine- 
qualities, |/c(z)| ^ 1 and |/c(z)| > 1 hold, respectively. Since by Assumption 1.4, 
/c = 0 is not a root of P{z . k) for | z | > 1, we obtain 

LEMMA 3.2 ([3], Lemma 4.3). For |z | ^ 1, f/ie r + p roofs k{z) of the basic 
characteristic function P{z , k) split into r inner roots with 0 < | k{z) \ ^ 1, and p 
outer roots with |/c(z)| ^ 1. 

We can now quote our main preliminary stability criterion in [3]: 

THEOREM 3.1 ([3], Theorem 4.2). ^Approximation (1.6) is stable if and only 
if for every z with | z | ^ 1 and all corresponding inner roots k(z) , the boundary 
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characteristic function E{z , k) satisfies 

i?(z, «(z)) # 0. (3.1) 

If k{z) is an inner root for z, |z | & 1, then by definition. 

P{z, k{z)) = 0. 

Thus, combining this fact with (3.1), we immediately obtain, 

COROLLARY 3.1. Approximation (1.6) is stable if and only if for every z 
with I z I > 1 and all corresponding inner roots k{z), 

n(z,/c(z))= |P(z,/c(z))| + \R{z,k{z))\ ts 0. (3.2) 

We shall need yet another result which describes the behavior of the inner 
roots for z = 1 : 

LEMMA 3.3 ([3], Lemma 5.1). If z = 1, then /c = 1 is not 071 inner root of 
P{z, k). 

This brings us to 

Proof of Theorem 2.1. Take any z, |z | > 1, and let k(z) be a corresponding 
inner root. If z 1. then by Lemma 3.2, we have 0 < | k{z) \ < 1; so the 
hypothesis in (2.1) implies (3.2). If on the other hand, z = 1, then by Lemma 3.3, 
k{z) = 1 is excluded as an inner root; so by Lemma 3.2 we have 
0 < I k(z) I < 1, k(z) ^ 1, eind again (2.1) yields (3.2). By Corollary 3.1 therefore, 
approximation 1.6 is stable, and the proof is complete. ■ 

Proof of Theorem 2.2. (i) If either the basic scheme or the boundeiry condi- 
tions are dissipative, then by Assumption 1.3(i) and Definition 2.1, either 

P(z. «) 0. |z I 2: 1. |/c| =1, « 1, 

or 


R{z,k)a0, \k\=1,kA^I\ 
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and (2.2a) follows. 

(iia) Since the boundary scheme is consistent with (1.5), it is at least first 
order accurate; so, in particular, it satisfies the ordinary zero order accuracy 
condition 

J=-T ff=0 j = -T 

which we equivalently write as 

P(z = 1, /c = 1) = 0. (3.3) 

By assumption, the basic scheme is two-level, so the function P(z, /c = 1) is a 
first degree pol 3 momial in z~' which, by (3.3), has a unique root z~^ = 1. Thus, 

P{z. K= 1) ^0. |z I > 1, z 1, 

and we obtain (2.2b). 

(iib) In the 3-level case we still have (3.3), but now P{z, k = 1) is a second 
degree pol 3 momial in z“* with real coefficients. By (3.3), z~^ = 1 is one of the 
roots, so the other root must be real as well; thus 

P{z, c = 1) 0, |z I = 1, z ^ ±1. (3.4) 

In addition, since the basic scheme satisfies the von Neumann condition 
(Assumption 1.2(i)), then 

P(z, K= 1) 0, |z|>l. (3.5) 

By (3.4) and (3.5) therefore, 

n(z,K = l)7«0, |z|^l, z jt ±1, (3.6) 

which together with the assumption in (2.5) gives us (2.2b). 

The proof of (iic) is identical to that of (iia), except that P{z, ic) is replaced 
by R{z, k), where by the zero order accuracy of the boundary conditions, we 
have R{z = 1, K = 1) = 0 instead of (3.3). 
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For the proof of (iid) we still have (3.5); and since the boundary conditions 
are three-level we may replace (3.4) by 

B{z, K= 1) 0. \z\ = l,Z7i±l. (3.7) 

By (3.5) and (3.7) we obtain (3.6) again, which together with (2.5) }delds (iid). 

(iii) Since the function z’^*i?(z , /c) is a polynomial in z and ic, it can be writ- 
ten as 

z«^ii?(z, k) = (z -Zi)- {z - z„) i?o(z. ic) (3.8) 

where i?o(z, k) is a polynomial in k whose leading coefficient is 1 and whose 
other coefficients are rational functions of z. By Assumption 1.6, the boundary 
characteristic function 

R{z,k)= 2 c,-(z)/c^ 
j=o 

is not the zero polynomial in k for any (fixed) z vrith |z | ^ 1. Thus, the Zj in 
(3.8) must satisfy 

|zj- 1 <1, ;• = 1 d. 

It follows that the leading coefficient of z’^*i?(z, k) in (3.8) is uniformly bounded 
away from zero for all |z | > 1. Hence, the roots k = k{z) of i?o(z, k), and there- 
fore those of R{z, k), are continuous functions of z for |z | ^ 1. 

Since the boundary conditions are assumed to be boundedly solvable, then 
by Lemma 2.1(i) we have 

T-M = 5^ 0, 0 < |k| < 1. (3.9a) 

Moreover, by (1.6c) 

T_i(0) = co(-i) / 0. (3.9b) 


So, combining (3.9a, b) we obtain 
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T-i(k) / 0. 1/c| < 1. (3.10) 

Since by hypothesis, the boundary conditions also satisfy the von Neumann 
condition, then for | z | > 1, the function R{z , k) does not vanish on the unit cir- 
cle |/c| = 1. Hence, by the continuity established above, the number of roots 
K = ic(z) of i?(z, k) satisf}dng |/c(z)| < 1, is independent of z, |z | > 1, and it 
thus equals the number of roots k, |/c| ^ 1, of the function 

k) = T-i{k). 

By (3.10) therefore, R{z, k) has no roots |/c| ^ 1 for |z | > 1, so we have 

R(z,k)t^0, |/c| < 1, |z|>l. (311) 

By (3.11), the roots k{z) of R{z, k) satisfy |k(z)| > 1 if |z | > 1. So for |z | > 1, 
these continuous roots satisfy |/c(z)| > 1, i.e., 

R{z , k) ^ 0, I K I < 1 , I z I ^ 1 , 

and (3.2c) follows. ■ 

In concluding the paper we make the following remarks: 

Remark 3.1. Our main results in [3] follow immediately from the present 
ones. Indeed, consider Theorem 3.3 of [3] where we assume that the basic 
scheme is dissipative, that the boundary conditions are boundedly solvable and 
satisfy the von Neumann condition, and that 

/?(z, K = 1) ^ 0, |z I = 1, z ¥■ 1. (3.12) 

With these assumptions. Theorem 2.2(i)(iii) gives (2.2a, c). Moreover, by 
Assumption 1.2(i) we have (3.5), which together with (3.12) yields (2.2b). Thus, 
Theorem 2.1' implies stability, and Theorem 3.3 of [3] follows. 

In Theorem 3.4 of [3] it is assumed that the basic scheme is arbitrary, that 
the boundary conditions are boundedly solvable and dissipative, and that (3.12) 
holds. Again, Theorem 2.2(i)(iii) )delds (2.2a, c), and (3.12) together with (3.5) 
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gives (2.2b): so Theorem 2.1’ implies Theorem 3.4 of [3]. 

Theorems 3.1 and 3.2 in [3] follow from Theorem 2.1’ and 2.2 with similar 
ease. 

Remark 3.2. Examples 2. 1-2.6 contain all the main examples in [3]. Exam- 
ple 2.3 (with a nondissipative basic scheme), and Examples 2.7-2.10, are new, 
however, in the sense that they could not have been handled by our old results 
in [3]. For examples 2.3, 2.7, 2.9 and 2.10 this is true since Theorems 3. 1-3.4 in 
[3] all require that either the basic scheme or the boundary conditions be dissi- 
pative. For Example 2.8, where the basic scheme is not necesseirily dissipative 
but the boundary conditions are, one may hope that Theorem 3.4 of [3] would 
help. But R{z = -1, /c = 1) = O; so (3.7) in [3] is violated, and Theorem 3.4 of [3] 
fails. 

Remark 3.3. Not all interesting cases are covered by our results. For 
instance, Gustafsson et al. showed in Theorem 6.1 of [5] that the Leap-Frog 
scheme in (2.8) together with the linear oblique boundary extrapolation 

+ Af) = Zvi{t) - vz{t - AO. 

provides a stable approximation. ¥6 easily find, however, that 

n(z =-l,K=-l)= \P{Z =-l, K=-l)| + \R{Z =-l.K=-l)| =0, 

so (2.1) may not hold. Consequently, Theorem 2.1 fails for this example, showing 
that our criteria are sufficient but not necessary for stability. 

Remark 3.4. Let us show that Assumption 1.5 is necessary for the stability 
of approximation (1.6). Indeed, if Assumption 1.5 fails, then at some z = zq with 
|zol — 1. all the Cj{z) in (1.12) vanish simultaneously. Hence, 

R{z = zo. 0=2 Cy(zo)/cJ = 0 ; (3.13) 

j=o ' 

that is, R{z = zq.k) is the zero polynomial in k. By Lemma 3.2, for z = zq, the 
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basic characteristic function P{z, ic) has r inner roots ic{zo) with 
0 < |/c(zo) 1 - !• I^or each of these roots (3.13) implies 

P{z = Zq, ic = k{zq)) = 0, 
and Theorem 3. 1 assures instability. 

For example, consider any dissipative basic scheme with boundary condi- 
tions determined by the two-level, zero order vertical extrapolation: 

v^{t + iit) = V = 0 r - 1. (314) 

Clearly, these explicit boundary conditions are boundedly solvable (Lemma 2.1), 
and they satisfy the von Neumann condition. Thus the hypotheses of Theorem 
2.2 hold, and Theorem 2.1' should have implied stability. As mentioned in Sec- 
tion 4 of [14], however, the approximation is unstable. The reason for this insta- 
bility is that the boundary conditions in (3.14) fail to satisfy Assumption 1.5, 
since 

2 icj(z)|z=i = |1 - = 0- 

J=0 

This example, though found unstable in the sense of Definition 3.3 of [5], is 
stable in the Iz sense described in [14]. 

Remark 3.5. It is known that the accuracy of the boundary conditions plays 
a decisive role in determining the overall rate of convergence of the difference 
approximation. In peirticular, Gustafsson [4] has shown that if approximation 
(1.6) is stable and the basic scheme is accurate of order 5, then an overall 6 rate 
of convergence is assured, if the boundary conditions are of 5 — 1 order of accu- 
racy at least. 

If we use horizontal or oblique boundary extrapolation as in (2.7) or (2.9), 
then obviously there is no restriction on accuracy at the boundary. In general, 
however, the question of maximal possible accuracy for the boundary conditions 
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can be quite difficult if the boundary conditions are to satisfy certain properties. 
For instance, it is not hard to see that the various boundary conditions in the 
examples of Section 2 (with the exception of horizontal extrapolation), are all 
stable in the sense of the definition in Assumption 1,2, with P{z, k) replaced by 
the characteristic boundary function R{z, k). Thus, it seems natural to raise 
the question of maximal accuracy for stable boundary conditions of type (1.6c). 
This question was answered by Iserles and Strang [7, Theorems 1 and 2] who 
treated the two-level case and showed that stable explicit boundary conditions 
of type (1.6c) may not be more than second order accurate, while implicit condi- 
tions of this t}q)e cem be constructed to arbitrary order of accuracy. 
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